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ABSTRACT 

The merger of two supermassive black holes is expected to produce a gravitational-wave sig- 
nal detectable by the satellite LISA. The rate of supermassive-black-hole mergers is intimately 
connected to the halo merger rate, and the extended Press-Schechter formalism is often em- 
ployed when calculating the rate at which these events will be observed by LISA. This merger 
theory is flawed and provides two rates for the merging of the same pair of haloes. We show 
that the two predictions for the LISA supermassive-black-hole-merger event rate from ex- 
tended Press-Schechter merger theory are nearly equal because mergers between haloes of 
similar masses dominate the event rate. An alternative merger rate may be obtained by invert- 
ing the Smoluchowski coagulation equation to find the merger rate that preserves the Press- 
Schechter halo abundance, but these rates are only available for power-law power spectra. We 
compare the LISA event rates derived from the extended Press-Schechter merger formalism 
to those derived from the merger rates obtained from the coagulation equation and find that 
the extended Press-Schechter LISA event rates are thirty percent higher for a power spec- 
trum spectral index that approximates the full ACDM result of the extended Press-Schechter 
theory. 
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1 INTRODUCTION 

Structure formation proceeds hierarchically, with small over- 
dense regions collapsing to form the first dark-matter haloes. 
These haloes then merge to form larger bound objects. The 
extended Press-Schechter (EPS) formalism provides a descrip- 
tion of "bottom-up" structure f ormation by combining t he Press- 
Schechter halo mass function jPress & Schec htei] l l974) with the 
halo merger rates derived by iLacev & Cold il993h . Since its in- 
ception, the EPS theory has been an invaluable tool and has been 
applied to a wide variety of topics in structure formation (see 
iBenson. Kamionkowski. & Hassarul2005l and references therein). 
Unfortunately, the Lacey-Cole merger-rate formula, which is 
the c ornerstone of EPS m erger theory, is mathematically inconsis- 
tent tee nson et all200a) . It is possible to obtain two equally valid 
merger rates for the same pair of haloes from the EPS formalism. 
These two merger rates are nearly equal when the masses of the 
two haloes differ by less than a factor of one hundred, but they di- 
verge rapidly for mergers between haloes with larger mass ratios. 
Consequently, any application of EPS merger theory gives two an- 
swers, and if the calculation involves mergers between haloes of un- 
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equal masses, the discrepancy between these two predictions may 
be large. 

Motivated by the ambiguity in the L acey-Cole merger rate, 
iBenson. Kamionkowski. & Hassanil i2005l hereafter BKH), pro- 
posed a method to obtain self-consistent halo merger rates. Since 
haloes are created and destroyed through mergers, the halo merger 
rate determines the rate of change of the number density of haloes 
of a given mass. By inverting the Smoluchowski coagulation equa- 
tion lSmoluchowskJ.191&) . .BKH find merger rates that predict the 
same halo population evolution as the time derivative of the Press- 
Schechter mass function. In addition to eliminating the flaw that 
resulted in the double- valued rates in EPS theory, the BKH merger 
rates by definition preserve the Press-Schechter halo mass distri- 
bution when used to evolve a population of haloes. The Lacey- 
Cole merger rate fails this consistency test as well, and the use of 
EPS merger trees has been constrained by this inconsistency (e.g. 
Menou. Haiman. & Narayanan 2001) . 

There are three limitations to the BKH merger rates. First, 
they are not uniquely determined because the Smoluchowski equa- 
tion does not provide sufficient constraints on the merger rate. The 
BKH merger rate is the smoothest, non-negative function that satis- 
fies the coagulation equation; it exemplifies the properties of a self- 
consistent merger theory, but it is not a definitive result. Second, 
the inversion of the Smoluchowski equation is numerically chal- 
lenging and solutions have been obtained only for power-law den- 
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sity power spectra. Finally, the BKH merger rates are derived from 
the Press-Schechter halo mass function rather than the mass func- 
tions obtained from N-body simulations I Sheth & Tormen \9^ % 
Ijenkins et alj200ll) . 

In this paper, we explore the possible quantitative conse- 
quences of our limited understanding of merger rates for one of the 
astrophysical applications of merger theory: the merger rate of su- 
permassive black holes. Since supermassive black holes (SMBHs) 
are believed to lie in the centre of all dark-matter haloes above 
some critical mass, halo mergers and SMBH mergers are inti- 
mately related. By considering only halo mergers that would re- 
sult in a SMBH merger, the EPS merger rates have been used to 
obtain SMBH merger rates taaehnelt 1994'; 'Menou et al."200l'; 
[Wvith e & Loeb 2003a,; .Sesana et al. 2004, 2005; Rhook & Wvithe 
[2005h . 

SMBH mergers are of great interest because they produce a 
gravitational-wave signal that may be detectable by the Laser Inter- 
ferometry Space Antenna (LISA), which is scheduled for launch 
in the upcoming decade. Consequently, EPS merger theory has 
been u sed to obtain es t imates for the SMBH merger event rate for 
LISA teaehnelt" 1994; "Menou et al."200l'; Wvithe & Loeb 2003a; 
fSesana et al. 2004, 2005; Rhook & Wvithe 2005) . In addition to 
their intrinsic interest as a probe of general relativity, there is hope 
that LISA'S observations of SMBH mergers would provide a new 
window into astrophysics at high redshifts. Wvithe & Loeb ( 2003a) 
used EPS merger theory to derive a redshift-dependent mass func- 
tion for haloes containing supermassive black holes and then used 
EPS merger theory to predict the LISA event rate that arises from 
this SMBH population. Since SMBH formation becomes more dif- 
ficult after reionization due to the limitations on cooling imposed 
by a hot intergalactic medium, the Wyithe-Loeb SMBH mass func- 
tion and corresponding LISA event rate are highly sensitive to the 
redshift of reionization. iMenou et alJfeOOll) used EPS merger trees 
to demonstrate that LISA observes more SMBH merger events 
when SMBHs at redshift 2; = 5 are only found in the most mas- 
sive haloes as opposed to being randomly distributed among haloes. 
iKoushiappas & Zentnen i2006) also used EPS merger trees to show 
that higher-mass seed black holes (Mbh ~ 10^ M© as opposed to 
Mbh ~ 10^ Mq) at high redshifts result in significantly higher 
LISA SMBH-merger event rates. Unfortunately, these ambitions 
of using LISA SMBH-merger event rates to learn about reioniza- 
tion and SMBH formation rest on the shaky foundation of extended 
Press-Schechter merger theory. 

We first review how the rate of mergers per comoving volume 
translates to an observed event rate in a ACDM universe and how 
the mass of the halo is related to the mass of the SMBH at its centre 
in Sections|2|and|3| In Section|4| we use the EPS formalism to de- 
rive an event rate for LISA. Throughout the calculation, we present 
the results derived from both versions of the Lacey-Cole merger 
rate. In Section|5] we explore the alternative merger-rate formalism 
proposed by Benson. Kamionkowski, & Hassani (2005). Since the 
BKH merger rates are only available for power-law density power 
spectra, it is not possible to use them to make a new prediction of 
the SMBH merger rate and the corresponding event rate for LISA. 
Instead, in Section|6| we use the event rates for power-law power 
spectra derived from the EPS and BKH merger theories to gauge 
how the LISA event rates may be affected by switching merger for- 
malisms. Finally, in Section^ we summarize our results and dis- 
cuss how these ambiguities in halo merger theory limit our ability 
to learn about reionization and supermassive-black-hole formation 
from LISA'S observations. 



2 COSMOLOGICAL EVENT RATES 

The merger of two supermassive black holes will produce a 
gravitational-wave burst. The observed burst event rate depends on 
the number density and frequency of black-hole mergers: the num- 
ber of observed gravitational-wave bursts per unit time (B) that 
originate from a shell of comoving radius R{z) and width AR is 



dB = (1 + zy^M{z) A-kR^ AR, 



(1) 



where N{z) is the SMBH merger rate per comoving volume as a 
function of redshift. The factor of (1 + 2:)^^ in equation results 
from cosmological time dilation. In equation Q, and throughout 
this article, we assume a flat ACDM universe. Given the relation be- 
tween comoving distance and redshift, AR — [c/H{z)] Az, equa- 
tion may be converted to a differential event rate per redshift 
interval. 



AB 

1z 



(1 



4n[R{z)fj^{z)c 
Hoy/nM{l + z)^ + nA 



(2) 



where SIm and Qa are the matter and dark-energy densities today 
in units of the critical density. 

The observed gravitational-wave burst rate from SMBH merg- 
ers is obtained by integrating equation J2j over the redshifts from 
which the bursts are detectable. LISA will be able to detect nearly 
all mergers of two black holes with masses greater t han 10'* Mp, 
and less than 10** Mq up to 2 ^ 9 taaehnell|[T994l ISesana et alJ 
[2005; Rhook & Wvithe 2003). Since more massive binary-black- 
hole systems emit gravitational radiation at lower frequencies and 
the observed frequency decreases with redshift, very distant [z ~ 
9) mergers of SMBHs with masses greater than 10^ Mp pro- 
duce signals below EISA's frequency window JSesana et all2005t 
Rhook & Wvithe 2005). However, the number density of 10* Mq 
haloes is exponentially suppressed at redshifts greater than four, so 
it is extremely unlikely that two black holes larger than 10** M0 
will merge at redshifts 2 Jl 4. Thus, the upper bounds on the rele- 
vant redshift and SMBH mass intervals are determined by the pop- 
ulation of supermassive black holes and not EISA's sensitivity. 



3 THE RELATIONSHIP BETWEEN HALO MASS AND 
BLACK HOLE MASS 

The transition from the rate of halo mergers to the rate of detectable 
SMBH mergers [Af{z) as defined in equation 0] requires a rela- 
tionship between the mass of a halo and the mass of the SMBH at its 
centre. Since LISA is sensitive to SMBH mergers at high redshifts, 
this A/bh — Mhaio relation must be applicable to high redshifts as 
well. 

Observations of galaxies out to 2 ~ 3 reveal a redshift- 
independent correlation between the m ass of the central black hole 
and the bulge velocity dispersion (Jc ^Ferrarese & Merrita 120001 
Gebhardt et al. 2000; Tremaine et al. 2002). The connection be- 
tween (7c and halo mass is mediated by the circular veloci ty v^. 
Using a sample of thirteen spiral galaxies. (Ferraresd i2002ah mea- 
sured a relationship between Vc and a^. Combining these measure- 
ments with the compiled relationship between SMBH mass and a^ 
presented by Ferrarese (2002b) reveals that observations are con- 
siste nt with a redshif t -indepe ndent A/bh oc v^ relation. 

IWvithe & Loebl y003b) proposed a mechanism for black- 
hole-mass regulation that would result in a A^bh oc v^ relation 
between central-black-hole mass and disc circular velocity for all 
redshifts. They postulated that a black hole ceases to accrete when 
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Figure 1. The masses of haloes that contain supermassive black holes of 
mass 10^, 10*, 10^ iqS and lO'^ Mq, according to the A/aH-A/halo re- 
lation proposed bylWyithe & Loeb ( 2003b) for a flat ACDM universe with 
f^M = 0.27. This relation is normalized to fit local observations and as- 
sumes that the disc circular velocity equals the virial velocity. 



the power radiated by the accretion exceeds the binding energy of 
the host galactic disc divided by the dynamical time of the disc. As- 
suming that the accretion disc shines at its Eddington luminosity, 
the black hole stops growing when 



Mbh = 1.9 X lO'' 



\0.07j V • 



Mr. 



(3) 



. 350 km s 

where Fq is the fraction of the radiated power which is transferred 
to gas in the disc. Setting Fq to 0.07 brin gs equation El in to agree- 
ment with the observations presented by Ferraresd i2002ah . 

The final step in the determination of a halo-black-hole-mass 
relation is to connect the circular velocity to the halo mass via the 
virial velocity I Bar kana & Loeb 2001). The simplest possible as- 
sumption is that the circular velocity of the disc equals the virial 
velocity of the halo. This assumption is made bv Wvithe & Loebl 
i2003d) . and we assume that Vc — «vir throughout this paper. 
However, different relations between Vc and Hvir have been pro- 
pos ed and can signifi cantly impact the final A/BH-.Afhaio relation 
('see lFerraresel2002a) . 

Assuming that Vc = tivir, the halo mass then becomes a 
redshift-dependent function of the mass of the central black hole: 



Mh 



1012 M, 



= 10.5 
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where SIm (2) the matter density in units of the critical density as a 
function of redshift, f2^ = Q,m{z — 0), and Ac is the nonlinear 
over-density at virialization for a spherical top-hat perturbation in 
a ACDM universe: 

Ac = IStt^ + 82[Qm{z) - 1] - 39[Qm{z) ~ if. (5) 

Figure \\\ shows the masses of haloes that contain supermassive 
black holes of several masses. Citing the fact that the largest 
haloes observed at low redshifts appear to contain galaxy clus- 
ters with no central black holes. Wyithe & Loeb (2003b) argue that 



supermassive-black-hole growth was complete by z ^ 1 and that 
local SMBH masses reflect the limiting values at that redshift. Con- 
sequently, when determining the mass of a halo that contains a 
black hole of a given mass, we use the z — 1 value of equation 
J4j for all redshifts less than one. 

Some calculations of the LISA SMBH-merger event rate im- 
pose a minimum halo virial temperature instead of a minimum 
black-hole mass when calculating the lowe r mass bound on haloes 
that contribute to the S MBH merger rate iWvithe & Loebll2003 J: 
iRhook & Wvithel2005^ . This constraint reflects the fact that super- 
massive black holes only form when the gas within dark-matter 
haloes can cool. However, the relation between virial temperature 
and virial mass ( Barkana & Loeb 2001) may be be used to elim- 
inate the halo mass in equation J3 in favour of the virial tem- 
perature. The redshift-dependent terms cancel, leaving a redshift- 
independent relation between black-hole mass and halo virial tem- 
perature: 



Mbh = (267 Mj 



-5/3 



iu 



Tvir 



X 104 K 



5/2 



(6) 



Therefore, defining the minimum halo mass by a minimum halo 
virial temperature is nearly equivalent to defining it by a minimum 
black-hole mass via equation J4j. For example, requiring that the 
halo's virial temperature be significantly higher t han the temper- 
ature of the intergalactic medium, Tvir <^ 10^ K iWvithe & Loehl 
;2003a), corresponds to imposing a minimum black-hole mass of 
2.6 X 10* M0. The only discrepancy occurs when z < 1, because 
we assume that the MsH-Mhaio relation is fixed for redshifts less 
than one, while Tvir is still redshift dependent. However, we shall 
see that nearly all SMBH mergers occur at redshifts greater than 
one, so this difference is negligible. 



4 EPS MERGER THEORY AND LISA EVENT RATES 

4.1 Review of EPS merger theory 

The first pillar of extended Press-Schech ter (EPS) merger theory i s 
the Press-Schechter halo mass function I'Press & SchechteJl974) . 
which gives the number of haloes with masses between M and 
M + dAI per comoving volume: 



duh 



alo 



dlnM 



2po 
ttA/ 



/diner \ 
V dlnM m) 



mJ a{M,z) 



exp 



2a^(M, 



,(7) 



where po is the background matter density today, ScoU is the criti- 
cal linear over-density for collapse in the spherical-collapse model, 
and (t(M, z) is the root variance of the linear density field at red- 
shift z in spheres containing mass M on average. In a ACDM 
universe, 5coii deviates slightly from its Einstein-de Sitter value 
of ~ 1.686 when the cosmological constant begins to domi- 
nate the energy density of the Un iverse (Kitavama & Suto 199g; 
' '^' IWeinbera & Kamionkowski 2003). In this work, the fitting func- 
tion obtained bv Kitavama & Suto (.1996) was used to approximate 
^coii- When calculating a{M, z), we assumed a scale-invariant pri- 
mordial power spectrum and we used the transfer function provided 
by Eisenstein & Hu ( 1998). 

The second pillar of EPS merger theory is the merger prob- 
ability function derived bv lLacev & Cold il993ft . which gives the 
probability that a halo of mass Mi will become a halo of mass 
A/f = All + M2 per unit time, per unit acquired mass: 



d> 
didMa 
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' a{Mi,z) 



1- 



-3/2 



X exp 



aHMi,z) 

1 1 



(8) 



In this expression, D{t) is the linear growth function, and a dot 
denotes differentiation with respect to time. 

Equation (SJ is usually interpreted as the differential probabil- 
ity that a given halo of mass Mi will merge with a halo of mass 
between AI2 and A/2 + dAf2 per unit time, per increment mass 
change. Thus equation J8) already includes information about the 
abundance of haloes of mass M 2 , but not the abundance of haloes 
of mass Ml. Following BKH* it is revealing to examine a differ- 
ent quantity, which does not differentiate between the two merging 
haloes: the rate of mergers between haloes of masses Mi and M2 
per comoving volume. 

Number of A/i + A/2 Mergers 



/?(A/i,A/2,t) 



di d (Comoving Volume) 
dn{Mi;t)\ ( d^p 



dA/i 



df dA'/2 



dA/idA/2.(9) 



The EPS self-inconsistency documented by iBKHl manifests itself 
here. Although R{M\, A/2, i) must be symmetric in its mass argu- 
ments by definition, equation J9j is not symmetric under exchange 
of A/i and A/2. 

The mass asymmetry of EPS merger theory becomes most 
transparent when one defines a new function: the merger kernel. 
From its definition, it is apparent that /?(A/i, A/2, t) should be 
proportional to the number densities of both haloes involved in 
the merger. Extracting this dependence defines the merger kernel 
Q(A/i,A/2,t): 



7?(A/i,A/2,t) 



/dn(A/i;t) 



/ dw(A/2 



t) 



dA/i J \ dA/2 
xQ(A/i,A'/2,t)dA/idA/2. 



(10) 



In addition to isolating the source of the mass-asymmetry in EPS 
merger theory, the merger kernel enters into the coagulation equa- 
tion which is inverted to obtain BKH merger rates, as described in 
Section|5| 

The EPS merger kernel Q[M\,M2) is the probability func- 
tion given by equation J8j divided by the number density of haloes 
of mass A/2 given by equation Q. In effect, EPS merger theory 
includes two distinct merger kernels, depending on the order of the 
mass arguments. Thus, we define two mass-symmetric merger ker- 
nels: (3m(A/i, A/2) equals the EPS merger kernel with the more 
massive halo as the first argument, while Ql(A/i, A/2) equals the 
EPS merger kernel with the less massive halo as the first argument. 
Figure |2| illustrates the differences in the merger kernels Qm and 
Ql. Note that neither QM(Afi, A/2) nor Ql(A/i, A/2) are viable 
candidates for the true halo merger kernel because they are not 
smooth functions of halo mass. They are useful because they ex- 
pose the ambiguities hidden in applications of EPS merger theory. 

In order to avoid double counting mergers when calculat- 
ing a merger rate, it is common to restrict one mass argument 
to be larger than the other. Using the standard expression for 
the Lacey-Cole merger probability function, as given by equation 
Jsj, in such calculations is equivalen t to us ing (5m(A/i,A/2) or 
Ql(A/i,A/2). Specifically, |Haehnei3 I1994I) effectively used Ql 

to predict an event rate for LISA, whilelWvithe & Loeb ( 2003a); 

I' 1 1 I "^ 

Rhook & W vithe (2005) effectively used Qm. Using the other ver- 
sion of the EPS merger kernel in either of these calculations would 
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Figure 2. The two EPS merger kemefs for z = 0. Here, Qm is the Lacey- 
Cole merger kernel with the more massive halo as the first argument, and 
Ql is the same kernel with the less massive halo as the first argument. 
Results are shown for a flat ACDM universe with Cm = 0.27, h = 0.72 
and (Tg = 0.9. 



have yielded different results, as we show in Section|42| More gen- 
erally, any application of the Lacey-Cole merger probability func- 
tion uses some mixture of Qm and Ql, and changing the mixture 
will change the result of the calculation. 

4.2 LISA event rates from EPS theory 

The rate of SMBH mergers per comoving volume follows from the 
rate of halo mergers per comoving volume given in equation J9j: 



N{z 



/ dn(Ml,z) ^ 
V dM / 



\ dM , 



xQ{Mi,M2,z)dMidM2, (11) 



where A/min is the minimum halo mass that contains a black hole 
massive enough to be detected when it mergers with a black hole 
of equal or greater mass. The factor of 1/2 acc ounts for the dou- 
ble counting of mergers. Some calculations fe.g. lRhook & Wvithd 
|2005) only include mergers between haloes with mass ratios less 
than three and so integrate A/2 from A/i/3 to 3A/i. This restric- 
tion is motivated by dynamical-friction calculations that indicate 
that when a halo merges with a halo less than a third of its size, 
it takes longer than a Hubble tiine for their central black holes to 
merge iCoIdI. Maver, & Governato 1999). However, recent numer- 
ical simulations indicate that this restriction may be too strict; when 
gas dynamics are included, SMBHs with host-galaxy-mas s ratios 
greate r than three merge within a Hubble time 1 Kazantzid is et alJ 
I2OO5I) . We do not impose any restrictions on the halo-mass ratios, 
so our event rates are upper bounds arising from the assumption 
that every halo merger in which both haloes contain a SMBH re- 
sults in a SMBH merger. 

Since LISA should observe mergers between two SMBHs 
with masses greater than lO** Mq out to redshifts of at least eight, 
we generally use this minimum black hole mass to determine 
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Figure 3. The rate of SMBH mergers per comoving volume where both 
merging black holes have a mass greater than 10"^ Mq, 10* Mq and 
10^ Mq. The sohd (dashed) hnes show the results when the first argument 
of the Lacey-Cole merger kernel is the more (less) massive halo. Results 
are shown for a flat ACDM universe with Q,yi = 0.27, h = 0.72, and 
0-8 = 0.9. 



Mmin- The corresponding rates of SMBH mergers per comoving 
volume are shown in Figure |3| as well as the rates which corre- 
spond to different choices for the minimum mass of a SMBH. Both 
versions of M are shown to illustrate the difference between the 
two Lacey-Cole merger kernels. The crimp in N{z) al z = 1 re- 
flects the transition from a constant A/min (evaluated at 2: = 1) to 
the redshift-dependent form given by equation J4j. 

Once the rate N[z) of SMBH mergers per volume is known, 
equation J2} may be integrated over redshift to obtain an event rate 
for LISA. Figure |4] shows the LISA event rate as a function of the 
minimum halo mass that contains a black hole large enough to emit 
an observable signal. Figure|5|shows the event rate as a function of 
the maximum redshift of a detectable merger. In Figure|5| A/min is 
the mass of a halo that contains a black hole more massive than 10^ , 
10*, or 10^ M0 as determined by the A/bh — Afhaio relation given 
by equation J4j. These rates correspond to the M results depicted in 
Figure|3| Examination of these results reveals that increasing Zmax 
beyond Zmax = 6 has little effect on the event rate when A/min is 
greater than 10^ M0, as is the case when equation J4j is used to 
obtain the value of Afmin which corresponds to a minimum black- 
hole mass of 10* M© . The levelling of the event rate for Zmax ^ 6 
indicates that SMBH mergers are very rare at higher redshifts and 
that the event rate is dominated by mergers that occur at redshifts 
2: ^ 6. Therefore, the upper bound on LISA'S sensitivity to larger 
SMBH mergers at high redshifts will have little effect on the event 
rate. 

The event rates shown in Figures |4l and [S] differ signif- 
icantly from those ca lculated by Wvithe & Loeb (2003a) and 
iRhook & WvitheH2005h . Our eve nt rates are generally mu ch higher 
than the event rates reported bv IWvithe & Loea J2003al) because 
we do not exclude mergers between haloes with mass ratios greater 
than three from our SMBH merger rate. The event rates calcu- 
lated by Rhook & Wyithe l 2005) are even lower because they do 
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Figure 4. The gravitational-wave event rate from SMBH mergers as a func- 
tion of the minimum halo mass that contains a SMBH large enough to 
produce a detectable signal when it merges. Mergers at redshifts up to 
^max were included in this rate, and the five pairs of lines correspond 
to Zniax = 2, 4, 6, 8, 10. The solid (dashed) lines show the results when 
the first argument of the Lacey-Cole merger kernel is the more (less) mas- 
sive halo. Results are shown for a flat ACDM universe with Qm = 0.27, 
h = 0.72, and o-g = 0.9. 
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Figure 5. The gravitational-wave event rate from SMBH mergers as a func- 
tion of the maximum redshift of a detectable merger. Only mergers in which 
both black holes have masses greater than the given lower bound are in- 
cluded. The solid (dashed) lines show the results when the first argument 
of the Lacey-Cole merger kernel is the more (less) massive halo. Results 
are shown for a flat ACDM universe with f^M = 0.27, h = 0.72, and 
as = 0.9. 
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not assume that all haloes contain galaxies. The one case where 
our event rates are not substantially higher than those derived by 
LRhook & Wvithe 1 2005) is when the minimum black-hole mass is 
taken to be very high (Mbh ^ 10® Mq). In that case, the mini- 
mum halo mass is so high that nearly all mergers involve haloes 
of similar masses ( Mhaio ^ 10^^ M^), a nd the galaxy-occupation 
fraction derived by Rhook & Wvithe ( 2005) indicates that nearly all 
haloes of this size contain galaxies for redshifts greater than three, 
so our event rate of 12 per yea r is very similar to the re sult of the 
more sophisticated treatment of iRhook & Wvithd J2005l) . ^ 

Event rates obtained from both versions of the EPS merger 
kernel are shown in Figures|4|and|5| The differences between these 
results reveal the type of mergers that dominate the calculation. 
As shown in Figure|5] the event rate obtained from Qm is slightly 
higher than the rate obtained from Ql, indicating that mergers 
with halo mass ratios less than 10^ are dominating the sum (see 
Figure |2j. For a constant value of Mmin = 10® Mq, the differ- 
ence between the two versions decreases as the maximum redshift 
increases, as shown in Figure |4| This convergence indicates that 
the contribution from mergers between haloes of greatly unequal 
masses to the event rate dwindles as redshift increases. Since the 
lower bound on halo mass is constant with redshift, a decrease 
in unequal-mass mergers reflects a decrease in the population of 
larger haloes. The Press-Schechter mass function implies that the 
largest halo that is common at a given redshift is given by the func- 
tion M-t{z), which is defined by the relation a{AI,, z) = 5cou{z). 
When M > Mt,{z), the exponential term in equation Q domi- 
nates, and the number density of such halos is exponentially sup- 
pressed. Mt{z) decreases with redshift, reflecting the fact that at 
early times, massive halos had yet to form. Due to the exponen- 
tial decline in the number density of haloes greater than AI,{z), 
there is an effective upper bound to the integrals in equation <11> . 
which defines M[z). This upper bound on halo mass follows A/* 
and is less than 100 times greater than A/min at the redshifts which 
dominate the merger rate. 

The relevant mass range may be quantified by considering the 
ratio. 



C(U) 



iil. dA^.]-^ dA^(3t 



V dA/i , 



^ dA/2 , 



Q(A^i,A^2) 



^ /m_ dM^ /m_ d^^^ (^) (3^) Q(M.,M,) ' 



where the ^-dependence of all quantities has been suppressed. Us- 
ing the standard Lacey-Cole merger kernel when evaluating C{U) 
is equivalent to using the arithmetic mean of Qm and Ql. Fig- 
ure ||| shows the values of [/ for C = 0.9, 0.95 and 0.99. Since 
the dominant halo mass range is so narrow, it is possible to find a 
power-law power spectrum that accurately approximates the value 
of a{M) over this mass range, and both the exact and the approx- 
imate o-(A/) are shown in Figure|6| This approximation will allow 
us to apply BKH merger theory to the calculation of LISA event 
rates in Section|6| 



^ When we attempted to reproduce the differential event rates calculated by 
[Haehnelt 1 1994), we found that our rates are roughly a factor of two lower 
After extensive review and two independent calculations, we were unable 
to find any errors in our analysis. 
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Figure 6. The halo mass range that dominates the rate of SMBH mergers 
per comoving volume. The thi'ee curves mai'ked with percentages define the 
upper bounds of mass ranges that account for 90%, 95% and 99% of A/". 
Here, Af^jin is the mass of a halo that contains a SMBH of mass 10* Mq . 
The dotted curves are plots of log cr{M) with arbitrary normalizations. Re- 
sults are shown for a flat ACDM universe with f^M = 0.27 and erg = 0.9. 
The dashed lines are plots of logcr(Af) for a power-law power spectrum 
with n = —2.1 and erg = 0.9843, which is the best linear fit to log cr over 
the mass range between Afmin and the 99% curve for 2 < 5. 



5 BKH MERGER THEORY 

5.1 Solving the coagulation equation 

A merger kernel that preserves the Press-Schechter (PS) halo mass 
distribution must satis fy the Smoluchowkski coagulation equation 
JSmoluchow ski 1916). which simply states that the rate of change 
in the number of haloes of mass AI equals the rate of creation of 
such haloes through mergers of smaller haloes minus the rate haloes 
of mass M merge with other haloes. Adopting the shorthand n(A/) 
for the PS halo number density per interval mass and suppressing 
the redshift dependence of all terms, the coagulation equation is 



di 



n{Ad) 



n{M')n{M - M')Q{M' , M - M') AAd' 



n{A4) / n{Ad')Q(M,M')dM', 



(12) 



where Q{Adi, M2, z) is the desired merger kernel. The first term 
on the right-hand side is the rate of mergers per comoving volume 
that create a halo of mass AI. The second term is the rate of mergers 
involving a halo of mass AI per comoving volume — these mergers 
effectively destroy haloes of mass M. 

|BKH numerically invert the coagulation equation for Q for 
power-law density power spectra Pik) ex fc". When the den- 
sity power spectrum is a power law, (j{AI) is proportional to 
j\^-(3+")/6 sijjgg jjjg redshift-dependence of the PS mass func- 
tion enters via the ratio 5coii(2:)/cr(A/, z) = (A//A/*)<^+"''''', 
the 2-dependence of the PS mass function may be eliminated by 
expressing the masses in units of Af, (z). For a judicious choice 
of time variables, differentiating the PS mass function introduces 
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no ^-dependence, and the coagulation equation becomes redshift- 
invariant. Consequently, the coagulation equation only has to be in- 
verted once, for the resulting merger kernel Q(Mi/M, , M2/M*) 
is applicable to all redshifts. This simplification is only possible 
when the power spectrum is a power law. For more complicated 
spectra, the coagulation equation will have to be solved at multiple 
redshifts. 

When they numerically solve the coagulation equation on a 
discrete grid, |BKH require that the merger kernel be symmetric 
in its two mass arguments. However, this restriction is not suffi- 
cient to determine Q uniquely from the coagulation equation. On 
an N X N mass grid, the coagulation equation becomes A'^ equa- 
tions for the TV possible values of Ad. Meanwhile, the symmetric 
Q matrix on the grid, Qij = Q{Mi, M-j), has N(N -|- l)/2 in- 
dependent components. To break the degeneracy, 'BKH! impose a 
regularization condition. By minimizing the second derivatives of 
Q, they find the smoothest, non-negative kernel that solves the co- 
agulation equation. 

5.2 BKH merger rates for power-law power spectra 

In Section l42l we demonstrated that the rate of SMBH mergers 
per comoving volume is dominated by mergers between haloes 
in a very limited mass range, as shown in Figure |6| The a{M) 
curves in Figure |S| show that it is possible to accurately approx- 
imate o-{M) over the relevant mass range as originating from a 
power-law power spectrum. We consider a power-law fit for a{M) 
that extends over all masses that fall within the 99% mass range 
at any redshift less than five. The fit has a lower mass bound of 
5.44 X 10^ Mq, which is the value of Afmin at z = 5, and ex- 
tends to a mass of 4.26 x 10^* Mq. Over this range, a{M) is best 
fit by spectral index n = —2.1 normalized so that erg = 0.9843, 
as shown by the dashed lines in Figure|S| This n = —2.1 power- 
law approximation of cr(i\f ) is accurate to within 16% over this 
mass range. We chose to fit the mass range for 2 ^ 5 because the 
SMBH merger rate peaks at redshifts less than five when the mini- 
mum black-hole mass is greater than 10* M0, as shown in Figure 
|3| Also, when the mass range is lowered, the best-fitting spectral 
index decreases, and BKH merger rates have not been obtained for 
n < -2.2. 

The density power spectrum enters the EPS merger kernel 
only through a(M), so any power-law approximation that accu- 
rately models o{M) for K'h, M2, and M{ = Mi + M2 will accu- 
rately model the Lacey-Cole merger kernel Q{Mi , M2 ,z)- Unfor- 
tunately, the same is not necessarily true for the BKH merger ker- 
nels obtained by inverting the coagulation equation. Since the coag- 
ulation equation [equation <12H involves integrals over all masses 
and is solved for all masses on the grid, the solution Q{M\ , M2 , z) 
is dependent on g{M) over all masses and not just the arguments 
of the kernel. Therefore, while the power-law approximation ac- 
curately reflects the full ACDM result for EPS merger theory, the 
BKH merger rates obtained for the same power law may differ 
greatly from the merger rates that solve the coagulation equation 
for a ACDM universe. However, since the coagulation equation has 
not been solved for a ACDM power spectrum, we compare the EPS 
merger rates to the BKH merger rates for the same power law. This 
comparison demonstrates how the BKH merger rates differ from 
the EPS rates, but should not be considered a definitive description 
of merger rates in a ACDM universe. 

BKH merger kernels for a power-law power spectrum with 
n — —2.1 were obtained by inverting the coagulation equation 
on a 91 X 91 grid of logarithmically-spaced M/M, values rang- 




-6-4-2 2 4 
Log[mymJ 



Figure 7. The two EPS merger kernels and the BKH merger kernel for a 
n = —2.1 power-law power spectrum at z = 0. Here, Qm is the Lacey- 
Cole merger kernel with the more massive halo as the first argument, and 
Ql is the same kernel with the less massive halo as the first argument. 
Results are shown for SIm = 0.27. h = 0.72 and erg = 0.98. The low- 
mass cut-off of the curves arises from the M/M^ ^ 10~^ bound on the 
BKH merger kernel. 



ing from 10"^^ to 3000. For M/Af* values greater than 10"*, 
the merger kernel values are not dependent on grid resolution, 
which indicates that the kernel is a numerically robust solution of 
the discretized coagulation equation for masses above 10^** A/,. 
The A'/BH-A^haio relation [equation 0] implies that SMBHs with 
masses greater than 10^ Mq and redshifts less than ten reside in 
haloes with masses greater than 10* M©, while the z = Q value of 
Af, for the n = —2.1 power-law power spectrum is 6 x 10^^ Mq. 
Therefore, for all haloes that contain SMBHs capable of producing 
a gravitational wave signal detectable by LISA, A//Af* ^ 10^^, 
so the lower mass bound on reliable kernel values is of no concern. 

Unfortunately, the same is not true for the upper bound on 
Af /Af, . The upper bound on the halo masses which contribute to 
the SMBH merger rate H in EPS theory, shown in Figure |6| ex- 
tends to MjM^ ^ 10^ for 2^5. However, extending the mass 
grid to higher values of M /M, introduces numerical noise that pre- 
vents the kernels from converging as grid resolution is increased. 
Therefore, we must extrapolate the BKH kernel to higher masses. 
We bilinearly extrapolate the logarithm of the kernel with respect 
to the logarithms of its mass arguments. When used to extrapo- 
late from a grid with M/KU < 100, this recovers the kernel to 
within a factor of two. Moreover, ignoring mergers of haloes with 
Af/Af» > 3000 only slightly affects the gravitational-wave event 
rate calculated from the BKH merger rates: the event-rate reduction 
is less than 3%. Therefore, the errors introduced by our extrapola- 
tion of the BKH merger kernel are negligible. 

The differences between the BKH merger kernel and both ver- 
sions of the EPS merger kernel are illustrated by FigureQ The BKH 
merger kernel is less than both EPS kernels when the masses of 
the merging haloes are similar, and the difference increases as the 
haloes get smaller. For mergers between haloes with mass ratios 
greater than 10^, the BKH merger kernel is nearly equal to the EPS 
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Figure 8. The rate of SMBH mergers per comoving volume where both 
merging black holes have a mass greater than 10^ Mq, 10'* Mq, or 
10^ Mq. The dotted line shows the EPS merger kernel for a ACDM power 
spectrum with erg = 0.9. The dot-dashed curves are the results derived from 
EPS theory for a power-law approximation with n = —2.1 and erg = 0.98. 
The dashed curves are the BKH results for the same power law and normal- 
ization. These results all assume a flat ACDM universe with Qm = 0.27 
and h = 0.72. 



G 
> 



10= 



102 

u 
CO 
CD 

a 10' 



10° 



I ' I ' I ' I 



"n^T 



Mbh > 10= M^ 



Mbh > 10^ M^ 



Mbh > 10^ M, 



///'' 



! S 

— ! S 
i i 

- i i . 
(' / / 

"i 



■EPS for CDM {a^ = 0.9) 
-EPS for n = -2.1 {a^ = 0.98) 

-BKH for n=-2.1 

I'i-M I I , I I I , I I I I I I I I I 



10- 

0123456789 10 

Maximum z of Events 

Figure 9. The gravitational-wave event rate from SMBH mergers as a func- 
tion of the maximum redshift of a detectable merger. Only mergers in which 
both black holes have a mass greater than 10^ Mq, 10* Mq, or 10^ Mq 
are included. The dotted line shows the EPS merger kernel for a ACDM 
power spectrum with erg = 0.9. The dot-dashed curves are the results de- 
rived from EPS theory for a power-law approximation: n = —2.1 and 
(Tg = 0.98. The dashed curves are the BKH results for the same power 
law and normalization. These results all assume a flat ACDM universe with 
Qm = 0.27 and h = 0.72. 



kernel with the least-massive halo as the first argument (Ql) for all 
masses. Therefore, for an n = —2.1 power-law power-spectrum, 
Ql comes closer to solving the coagulation equation than Qm ■ 



6 COMPARISON OF LISA EVENT RATES FROM BKH 
AND EPS MERGER THEORIES 

Since the BKH merger kernels for haloes of nearly equal masses are 
smaller than the EPS kernels for the same spectral index, applying 
EPS merger theory may over-estimate the LISA event rate. Figure 
|8| shows the rate M of SMBH mergers per comoving volume for 
the power-law model discussed in the previous Section. For com- 
parison, the EPS results for a ACDM universe are also shown as 
dotted curves (these are the arithmetic means of the corresponding 
solid and dashed curves in Figure |3j. However, it is important to 
remember that although the power-law models may accurately ap- 
proximate the ACDM results in the EPS theory, the same should 
not be assumed for the BKH merger rates. The BKH merger rates 
should only be compared to the EPS rates for the same power law. 
The discrepancy between the power-law EPS results and the 
ACDM curves at high redshifts is attributable to the power-law halo 
number density, which is much greater than the ACDM halo num- 
ber density for masses below IQi^i^ Mq at these redshifts. The same 
mass function is used to calculate the merger rate in BKH theory, 
so when the power-law merger rate is higher than the ACDM rate 
in EPS theory, it is reasonable to assume that the same is true for 
the rate derived from BKH theory. Figure|8|also shows that the pre- 
dictions for the SMBH merger rate from the BKH and EPS merger 
theories diverge with increasing redshift. In Section l4!2l we showed 
that as redshift increases, nearly equal-mass halo mergers dominate 



the event rate. The differences between the BKH merger kernel and 
the EPS kernel are greatest when the masses of the merging haloes 
are nearly equal, so as these mergers dominate the event rate at high 
redshifts, the BKH and EPS merger rates diverge. 

Figure |9| illustrates the potential consequences BKH merger 
theory has for the SMBH merger event rate observed by LISA. The 
difference between the BKH and EPS merger kernels for the same 
spectral index leads to a fairly substantial difference in the resulting 
event rates for LISA. For realistic values of the maximum redshift 
of a detectable merger (zmax ^ 5), the EPS prediction is about 
thirty percent higher than the BKH prediction for the n = —2.1 
power-law approximation. If the BKH merger kernel for a full 
ACDM power spectrum preserves the ratio of the BKS and EPS 
event rates for this spectral index, the LISA event rate from SMBH 
mergers would be reduced as well. |Rhook & Wvithe ( 2005) used 
EPS merger theory to predict that LISA will have approximately 15 
SMBH-merger detections per year at a signal to noise greater than 
five (they only consider mergers with Mbh ^ 10^ Mq). These 
comparisons of EPS and BKH event rates indicate that EISA's event 
rate may be closer to ten, with all other assumptions held fixed. 



7 SUMMARY AND DISCUSSION 

The EPS merger theory used to predict supermassive-black-hole 
merger rates is mathematically inconsistent because it contains two 
merger rates for the same pair of haloes. When the EPS formal- 
ism is used to derive supermassive-black-hole merger rates and the 
corresponding event rate for LISA, there are two potential results; 
the EPS predictions are ambiguous. We have found that mergers 
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between haloes whose masses differ by less than a factor of 10^ 
dominate the SMBH merger rate, even when all mergers between 
SMBH-containing haloes are included. The difference between the 
EPS merger rates for mass ratios in this range is small, so the two 
merger rates predicted by EPS theory are nearly equal. 

The concordance between the two EPS predictions for the 
SMBH merger rate is an artifact of the relative paucity of haloes 
with masses larger than 10^^ Mq. It is not an indication that the 
EPS merger formalism may be trusted to give realistic merger 
rates. In addition to its mass-asymmetry, the Lacey-Cole merger 
rate fails to give the same evolution of the halo population as 
the time derivative of the Press-Schechter mass function. Both of 
these flaws justify the search for a new theory of halo mergers. 
[Benson. Kamionkowski. & Hassani (2005) (BKH) inverted the co- 
agulation equation to find merger rates that preserve the Press- 
Schechter halo mass function for power-law power spectra. They 
found that these merger rates differ significantly from the EPS rates 
for the same power spectrum. 

The limited range of halo masses that contribute to the SMBH 
merger rate makes it possible to find a power-law power spec- 
trum that accurately fits the mass variance a{M) in this region. 
We consider such a power-law approximation with spectral index 
n = —2.1. Since the EPS merger formula depends only on the val- 
ues of (t(A/) for the two halo masses that are merging and the mass 
of the resulting halo, the power-law approximation accurately de- 
scribes the result obtained from the ACDM power spectrum. The 
same correspondence cannot be assumed for the BKH merger rates 
because they are dependent on a{M) at all masses. 

Nevertheless, it is illuminating to compare the SMBH merger 
rates derived from BKH merger theory to those derived from EPS 
theory for the same spectral index. When n = —2.1, the BKH 
merger rates are lower than the corresponding EPS rates for nearly 
equal-mass halo mergers, which dominate the rate of SMBH merg- 
ers. This discrepancy is a clear demonstration of how the EPS rates 
fail to solve the coagulation equation and therefore fail to preserve 
the PS halo mass function. It also indicates how BKH theory may 
predict a different SMBH-merger event rate for LISA, since the 
difference in merger rates results in an equally large difference 
in event rates. Comparing the event rates derived from EPS and 
BKH merger theories for this spectral index indicates that the LISA 
event-rate predictions that employ EPS merger theory may over- 
estimate the event rate by thirty percent. 

Fortunately, the ambiguity carried into the SMBH-merger 
event-rate predictions for LISA from the uncertainty surrounding 
halo merger theory does not appear to immediately preclude ex- 
tracting information regarding reionization or blac k-hole forma- 
tion fro m LISA'S observations of SMBH mergers. Wvithe & Loebl 
J2003g) showed that the LISA SMBH-merger event rate with reion- 
ization occurring at 2; = 7 is about 2.4 times higher than if reion- 
ization occurred earlier, at z — 12. This difference is larger than the 
uncertainties in the event rate revealed by our comparisons of BKH 
and EPS predictions, so it may be possible to constrain the reion- 
ization redshift from the LISA SMBH-merger event rate without 
a definitive theory of halo mergers. The thirty-percent uncertainty 
implied by these halo-merger-theory comparisons is also less than 
the differe nce in event rates for different SMBH seeding found by 
iMenou et: al. (2001). However, a thirty-percent uncertainty in the 
SMBH-merger rate will significantly loosen the constraints EISA's 
observations of SMBH mergers could place on reionization and 
SMBH formation. More concerning is the fact that there is no guar- 
antee that the merger rate derived from the merger kernel that satis- 



fies the coagulation equation for a ACDM universe does not differ 
from the EPS merger rate by more than thirty percent. 

Clearly, solving the coagulation equation for a ACDM 
power spectrum is imperative. Any application of extended Press- 
Schechter merger theory to astrophysical phenomena has a flawed 
foundation and the resulting predictions are unreliable. Specifi- 
cally, we have shown that the differences between EPS merger the- 
ory and BKH merger theory for power-law power spectra indicate 
that switching merger theories could significantly alter the LISA 
SMBH-merger event rate. This theoretical uncertainty should be 
resolved before EISA's measurements of SMBH merger rates are 
used to constrain cosmological models. 

We thank lonathan Pritchard for useful discussions. ALE is 
supported by a NSF Graduate Fellowship. AIB is supported by 
a Royal Society University Research Fellowship. This work was 
supported at Caltech by DoE DE-FG03-92ER40701 and NASA 
NNG05GF69G. 



REFERENCES 

Barkana R., Loeb A., 2001, Physics Reports, 349, 125, [arXiv:astro- 

ph/00 10468] 
Benson A. J., Kamionkowski M., Hassani S. H., 2005, MNRAS, 357, 847, 

[BKH] [arXiv:astro-ph/0407136] 
Colpi M.. Mayer L., Governato R, 1999. ApJ, 525, 720, [arXiv:astro- 

ph/9907088] 
Eisenstein D. J., Hu W., 1998, ApJ, 496, 605, [ai-Xiv:astro-ph/97091 12] 
Ferrarese L., 2002a, ApJ, 578, 90, [ai-Xiv:astro-ph/0203469] 
— , 2002b, Current High Energy Emission Around Black Holes, Lee C, ed. 

Singapore: World Scientific, [arXiv:astro-ph/0203047] 
Ferrarese L., Merritt D., 2000, ApJ, 539, L9, [ai-Xiv:astro-ph/0006053] 
Gebhardt K., Richstone D., Kormendy J., Lauer T. R., Ajhar E. A., Bender 

R., Dressier A., Faber S. M.. Grillmair C, Magorrian J., Tremaine S., 

2000, AJ, 119, 1157, [arXiv:a.stro-ph/9912026] 
Haehnelt M. G., 1994, MNRAS, 269, 199, [arXiv:astro-ph/9405032] 
Jenkins A., Frenk C. S., White S. D. M., Colberg J. M., Cole S., Evrard 

A. E., Couchman H. M. R, Yoshida N., 2001, MNRAS, 321, 372, 

[arXiv:astro-ph/0005260] 
Kazantzidis S., Mayer L., Colpi M., Madau R, Debattista V. R, Wadsley 

J., Stadel J., Quinn T., Moore B., 2005, ApJ, 623, L67, [arXiv:astro- 

ph/0407407] 
Kitayama T., Suto Y., 1996, ApJ, 469, 480, [arXiv:astro-ph/9604141] 
Koushiappas S., Zentner A., 2006, ApJ, 639, 7, [arXiv:astro-ph/0503511] 
Lacey C, Cole S., 1993, MNRAS, 262, 627 
Menou K., Haiman Z., Narayanan V. K., 2001, ApJ, 558, 535, [arXiv:astro- 

ph/0101196] 
Press W. H., Schechter P, 1974, ApJ, 187, 425 
Rhook K. J., Wyithe J. S. B., 2005, MNRAS, 361, 1145, [arXiv:astro- 

ph/0503210] 
Sesana A., Haardt F., Madau P, Volonteri M., 2004, ApJ, 611, 623, 

[arXiv:astro-ph/0401543] 
— , 2005, ApJ, 623, 23, [ai-Xiv:astro-ph/0409255] 
Sheth R. K., Tormen G., 1999, MNRAS, 308, 119, [arXiv:astro- 

ph/9901122] 
Smoluchowski M., 1916, Physik. Zeit., 17, 557 
Tremaine S., Gebhardt K., Bender R., Bower G., Dressier A., Faber S. M., 

Filippenko A. V., Green R., Grillmair C, Ho L. C, Kormendy J., Lauer 

T. R.. Magorrian J., Pinkney J., Richstone D., 2002, ApJ, 574, 740, 

[arXiv:astro-ph/0203468] 
Weinberg N. N., Kamionkowski M., 2003, MNRAS, 341, 251, [arXiv:astro- 

ph/0210134] 
Wyithe J. S. B., Loeb A., 2003a, ApJ, 590, 691, [arXiv:astro-ph/0211556] 
— , 2003b, ApJ, 595, 614, [arXiv:astro-ph/0304156] 



© 0000 RAS, MNRAS 000, 000-000 



